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A model for nonequilibrium dynamical mean-field theory is constructed for the infinite dimensional 
Hubbard lattice. We impose nonequilibrium by expressing the physical orbital as a superposition of a 
left-(L) moving and right- (i?) moving electronic state with the respective chemical potential /hl and 
fiR. Using the second-order iterative perturbation theory we calculate the quasi-particle properties 
as a function of the chemical potential bias between the L and R movers, i.e. $ = /xl — I^r- The 
evolution of the nonequilibrium quasi-particle spectrum is mapped out as a function of the bias and 
temperature. The quasi-particle states with the renormalized Fermi energy scale eqp disappear 
Sit ^ ^ Sqp in the low temperature limit. The second-order perturbation theory predicts that 
in the vicinity of the Mott-insulator transition at the Coulomb parameter U — Uc^ there exist 
another critical Coulomb parameter Ud (< Uc) such that, for Ud < U < Uc, quasi-particle states 
are destroyed discontinuously when (sqp)'^ ~ a{7TkBTc)^ + b^c with the critical temperature Tc and 
the critical bias $c- 
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In recent years significant experimental progress has 
been made in the fabrication of sophisticated electronic 
heterostructures utilizing strongly correlated electronic 
materials. These systems have given rise to the discovery 
of rich novel phenomena. These include: ballistic trans- 
port of electrons through heterostructures of supercon- 
ductor^, ferromagnet^ magnetic tunneling junction^^ 
and oxide^^'^^'. With this vast array of strongly cor- 
related heterostructures it is very important to theoret- 
ically understand how the strongly correlated materials 
in the bulk limit will behave under nonequilibrium con- 
ditions. In this work we formulate the nonequilibrium 
problem in the lattice and present significant progress 
towards a more complete understanding of strongly cor- 
related lattices out of equilibrium. 

We are interested in heterostructures of strongly corre- 
lated materials under a finite source-drain bias where the 
voltage drop occurs mainly at the interface of the s trongly 
correlated material and the source / drain lead^l^. Inside 
the strongly correlated material the transport is ballistic 
and driven by the momentum distribution of imbalanced 
chemical potentials. Our focus here is to investigate the 
evolution of strongly correlated quasi-particle spectra due 
to nonequilibrium driven by an imbalance in the chemical 
potentials. In this work we gain qualitative understand- 
ing on how the enhanced dephasing by the new particle- 
hole decay channel due to nonequilibrium modifies the 
spectral properties. 

One of the most practical and powerful theoretical 
techniques for studying strongly correlated lattices is the 
dynamical mean-field theory (DMFT)^^. Within this 
approach the self-energy of the strongly correlated lat- 
tice becomes momentum-independent and as a result the 
problem is reduced to solving a self-consistent interact- 
ing impurity model. This inherent simplicity of DMFT 
makes it a very attractive tool, especially with the large 
quantity of exact and perturbative impurity solvers such 



as the Hirsch-Fye quantum Monte Carlo method , the 
numerical renormalization group^^, and the non-crossing 
appr oximat iorP^I^. 

As an extension of DMFT to nonequilibrium situ- 
ations, Okamoto has theoretically studied the nonlin- 
ear transport and spectral properties of metal-Mott 
insulator-metal heterostructures where a bias voltage 
is applied across the Mott insulatoi*^^ using layered- 
DMFT^^. There he combined the layered-DMFT tech- 
nique with the Keldysh Green function approach and in- 
voked the non-crossing approximation as the impurity 
solver. Another system recently studied is the DMFT 
limit of the Hubbard model with the nonequilibrium 
driven by a uniform high electric fielcf^. Others have 
looked at the DMFT limit of the Falicov- Kimball model 
in the presence of a uniform time-dependent electric field 
and examined the transient current, quenching of the 
Blo ch o scillations, and evolution of the spectral func- 
tiorPI^. 

In our analysis of the problem we choose to start 
from the metallic state and analyze how the Fermi liq- 
uid is renormalized and eventually destroyed with multi- 
ple chemical potentials. We find that the quasi-particles 
strongly depend upon the strength of the chemical po- 
tential difference ^. At zero $ it has long been known 
that the system undergoes a metal-insulator transition at 
U = Uc^^. The second order iterative perturbation the- 
ory (IPT) approximation gives Uc = Uc2 — 3.31) at zero 
temperature, where D is the half bandwidth. As the tem- 
perature is raised the transition remains distinct down 
to Uci — 2.6D. In nonequilibrium, the second order IPT 
shows that the quasi-particles are destroyed abruptly by 
the chemical potential bias in the region Ud < U < Uc, 
where Ud — 2.3D. When U < Ud the system exhibits 
a smooth crossover from a system of well-defined quasi- 
particle states to that without quasi-particle excitations. 

We start with a non-interacting (i-dimensional tight- 
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binding Bethe lattice. The tight-binding solutions are 
grouped evenly into left- (L) and right- {R) movers with 
chemical potentials jiiL = *^*/2 and jj^n = — <^>/2, respec- 
tively. Since we do not study transport properties in this 
work, it is not important how the physical site is decom- 
posed into L, R movers. The original site orbital dl^ at 
site i with spin a is written as a superposition of L and 



R orbitals (i 1 / ^ • according to 



dl 



)• 



(1) 



The interaction terms and any observable quantities 
ought to be given in the physical basis d]^. The nonequi- 
librium DMFT lattice is depicted in Figure [l] 
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FIG. 1: The Bethe lattice where the (physical) d-staie is given 
as a superposition of a left-(L) and a right- (i?) moving state 
with the corresponding chemical potentials /jl = ^/2 and 
flR = -$/2. 



of DMFT is applied within the physical basis where 
the nonequilibrium statistics are imposed on the non- 
interacting Green functions. 

This Hamiltonian is solved using the iterative pertur- 
bation theor};^^. The IPT has been extensively utilized 
and well-established in equilibrium DMFT and we an- 
ticipate that the second order IPT gives a qualitatively 
reasonable description of the quasi-particle destruction 
in nonequilibrium. However we caution that the second 
order perturbation theory may not capture the exact na- 
ture of the quasi-particles in nonequilibrium. For exam- 
ple, it is known that second-order perturbation theory 
does not properly describe the splitting of the Kondo 
peak in quantum dot devices with a finite voltage bias, 
as seen in the fourth-order perturbation theor}/^^ and in 
nonperturbative solvers^^^^. Therefore, it still remains 
to be seen through nonperturbative calculations which of 
the following results will hold up. 

Now we outline the DMFT self-consistent routine. The 
calculation is performed in real-time as opposed to real- 
frequency because close to the transition of a system with 
well defined quasi-particle states to a system with none 
the bath spectral function becomes sharply peaked at the 
Fermi energ}/^^, thus making the real- frequency calcula- 
tion much more difficult. A typical time bandwidth we 
use is ~ 3000. In the non-interacting limit the (even) 
bath green function is given by 



(4) 



The model with which we focus is the DMFT limit 
{d oo) of the Hubbard model H = Hq -\- V , 



Ho 



<i,j> a=L,R;a 



(2) 
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where /^(e) = /(e - |) and /i^(e) = /(e -f |) are the 
Fermi functions for the L and R movers. D{e) is the 
non-interacting density of states (DOS), which we have 
taken to be semi-circular D{e) = :^^\/D'^ — e^. In the 
following, D = 1 is used as the unit of energy. 

Within IPT the self energy is calculated to second- 
order in U. Using the Langreth theorem^^ we may cal- 
culate the lesser and greater self energies according to 



The non-interacting tight-binding part of the Hamilto- 
nian 7io takes into account the hopping of electrons be- 
tween nearest neighbor lattice sites of the Bethe lat- 
tice. The interaction term of the Hamiltonian V is 
the particle-hole symmetric on-site Coulomb interaction 
of strength U. Both the non-interacting and interact- 
ing terms of the Hamiltonian are homogeneous in real 
space, which is necessary for our application of single- 
site DMFT. The Hamiltonian may be further simpli- 
fied by transforming the L and R basis into the even 
(E) and odd (O) superposition of the basis. The cre- 
ation operators for the even and odd electron orbitals 
are d^^^Q.^ = {d]^ ± ^^)/V^. The interaction is 
only given by the even (physical) basis and the non- 
interacting part is decoupled into the even-odd parts 
as Ho = -tE<i,j> Ea^Ej-A^iM^Jcr + ^.c). There- 
fore, the time-evolution of even (physical) and odd states 
are completely decoupled and the perturbation theory 



U'[g^{t)rg^{-t). 



(5) 



The lesser and greater self energies are then Fourier trans- 
formed to real frequency. The lesser (greater) self energy 
gives us the particle (hole) spectral weight as a func- 
tion of frequency for the retarded self-energy. Using 
the Keldysh Green function relation — T^f^^{uj) = 

^i^t(^)~^i^t(^)' express the retarded self energy 

in the spectral form 



sr„t(^) 
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Im[S<,(e)-S>,(e)] 
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The nonequilibrium DMFT self consistent equations for 
the interacting and bath Green function are therefore 



de- 



D{e) 



(7) 
(8) 
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where po{uj) = — 7r~-^Im[^^(cj)] is the new bath DOS. For 
the new Keldysh bath Green functions we replace the 
non-interacting DOS in Eq. Q with po(^)- The mo- 
mentum independence of the self energy and the trans- 
lational invariance ensure that the wave vectors remain 
good quantum numbers. We iterate these equations until 
the bath Green function and interacting Green function 
converge. 
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FIG. 2: (color online) Interacting local spectral functions for 
U = 2, inverse temperature p = 300 and half bandwidth 
D = 1 plotted for a range of chemical potential biases, 
At this value of the Coulomb interaction U the quasi-particle 
peak is destroyed continuously as $ is increased. 

From previous works in the equilibrium systemF^, it 
is known that by increasing the temperature the quasi- 
particles are destroyed when the thermal fluctuations sur- 
pass a certain low-energy scale. This low-energy scale 
is given by the renormalized Fermi energy €qp which 
we define as the half- width at half maximum (HWHM) 
of the quasi-particle peak and is similar to = ZD 
where Z = [1 — dujRe[T^^ {u)]\^=o]~^ is the quasi-particle 
weight, as found in Georges et.afi^. As the temperature 
is increased at zero bias <l> = 0, the quasi-particles are 
destroyed when kpTc ~ Sgp where Tc is the critical tem- 
perature. £qp is the quasi-particle half bandwidth at 
zero bias. 

Figures [2][4] show the evolution of quasi-particle spec- 
tra as a function of bias <l>. The destruction of the Fermi 
liquid is easily understood by enhanced dephasing of par- 
ticles and holes, which is a consequence of the opening 
up of phase space available for the particles/holes to scat- 
ter into through <l> and by thermal fiuct nations. Both of 
these effects lead to a finite lifetime for the electron at 
the Fermi energy. To lowest order in cj, T and ^ the 
imaginary part of the self energy obeys the relation^^ 



lim Im[S^(a;) 



(9) 



The interacting spectral function p{ijo) is calculated both 
as a function of the Coulomb interaction U and the 
applied chemical potential bias ^. For U < Ud with 




FIG. 3: (color online) Interacting local spectral functions for 
U — 2.3, 2.5, 2.6, and 2.8 at inverse temperature [3 — 300 and 
band width D = 1. These plots depict the sudden disappear- 
ance of the quasi-particle peak at the critical voltage $c. The 
spectral functions are plotted at $ = 0, at $ just before the 
quasi-particle destruction, and at $ = $c- 



Ud 2.31), the destruction of the quasi-particle peak is 
continuous, and exhibits similar behavior to the metal- 
insulator transition in the crossover regime {U < Ud) of 
equilibrium DMFT. Our results at = 2D{< Ud) are 
shown in Figure [2] As the chemical potential bias is in- 
creased the quasi-particle weight smoothly shifts towards 
the upper and lower Hubbard bands, and the quasi- 
particle peak disappears when ^ ~ ^qp- 

For Ud < U < Uc the quasi-particle peak is discon- 
tinuously destroyed at a critical chemical potential bias 
^c- In Figure [3] we have plotted the interacting spectral 
functions for U/D = 2.3, 2.5, 2.6, and 2.8 to illustrate the 
sudden disappearance of the quasi-particle peak. The 
spectral function for each value of U is shown for <l> = 0, 
<l> just before the transition, and finally at the transition 
<l> = This discontinuous transition in nonequilibrium 
by bias is reminiscent of the discontinuous transition by 
temperature in the equilibrium DMFT. 

In Figure |4] the scaled quasi-particle energy £qp/sqp 
is plotted versus the scaled chemical potential bias ^/^c- 
For ^ < 0.4^c the curves for U/D = 2.3, 2.5, 2.6 and 2.8 
scale onto a single curve. At <l> = it is clearly visible 
that the quasi-particle peak disappears discontinuously 
for each value of U. 

In Figure |5] we take a closer look at the temperature 
dependence of the transition in the discontinuous region. 
At U = 2.8, which is well into this region, the quasi- 
particle disappearance remains discontinuous for the en- 
tire range of temperatures. Therefore the discontinuous 
nature of the transition is very robust and is not effected 
by the temperature. The inset yields the dependence 
of the critical bias upon temperature. Motivated by the 
form of the imaginary part of the self energy [Eq. ([9|], we 
anticipate that the destruction of the quasi-particles will 
occur when (sgp)'^ ~ aiTikpTc)'^ + at critical values 
for the temperature Tc and the bias In Figure [5] we 
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FIG. 4: (color online) Quasi-particle energy sqp as a function 
of the chemical potential bias O at /3 = 300 and D — 1. The 
quasi-particle energy is scaled to the quasi-particle energy at 
zero bias, £qp, and the chemical potential bias is scaled to the 
critical bias at which the quasi-particle peak is destroyed, Oc- 
For small bias, O < O.40c, the quasi-particle energies scale to 
a single curve. In this range of values for U the disappearance 
of the quasi-particle peak is strongly discontinuous. 
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FIG. 5: (color online) Quasi-particle energy £qp as a func- 
tion of the chemical potential bias O at [/ = 2.8 and D — 1 
for different values of the inverse temperature f3. As the tem- 
perature is increased, the disappearance of quasi-particles re- 
mains discontinuous and the critical bias Oc is lowered. The 
inset gives the temperature dependence of the critical bias. 
The data (black circles) was fit (blue curve) to the function, 
{^Qpf ^ a{7TkBTcf + b^l, where a = 1.1 and b = 0.54. 



fit versus Tc to this function and find that a = 1.1 
and b = 0.54. Considering the nonlinear effects in the 
numerical results, the obtained values a and b are in rea- 
sonable agreement with the estimates 1 and | based on 



The critical bias and renormalized Fermi energy 
Sqp are plotted versus U in Figure |6| at = 300. The 
onset of the discontinuous transition is marked dit Ud — 
2. 3D and the Mott-insulator transition at <l> = is given 
at Uc ^^3.3D. We see that the transition occurs at ~ 



0.8 - 




Eq. m 



FIG. 6: The critical bias Oc is plotted as a function of the 
Coulomb interaction U at inverse temperature P = 300 and 
bandwidth D = 1. We have also plotted the renormalized 
Fermi energy £gp at O = 0. In the region Ud < U < Uc 
the quasi-particle peak is destroyed discontinuously at Oc- 
When U < Ud the quasi-particle peak disappears smoothly 
with increasing chemical potential bias. At this temperature 

Oc - Sqp. 



Sqp which further justifies the function used for fitting 
<l>c versus ksTc of Figure [5] 

In this work we have presented a theoretical model 
for nonequilibrium dynamical mean-field theory and per- 
formed the iterative perturbation calculations within 
the dynamical mean-field theory. In a lattice where 
the statistics are determined by the presence of multi- 
ple chemical potentials, the quasi-particle properties are 
strongly dependent upon the strength of the Coulomb in- 
teraction and the chemical potential bias, ^. For U < Ud 
with Ud — 2. 3D for the bandwidth D^ the quasi-particle 
disappears smoothly with increasing chemical potential 
bias. The disappearance of the quasi-particles is caused 
by the opening up of phase space for the electrons to scat- 
ter into at the Fermi energy, resulting in a finite lifetime 
for the electrons. The perturbation theory also predicts 
that in the region Ud < U < Uc, where Uc — 3. 3D marks 
the Mott-insulator transition in equilibrium, the quasi- 
particle particles may be destroyed discontinuously at a 
critical bias, ^c- 
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